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We have performed similar studies for an alternate can- 
didate for the order parameter symmetry. As shown in 
Ref [15], when the c- axis modulation of the spin fluctua- 
tion spectra are included in the gap equation, the two in 
phase rf-states of the bi-layer unit, compete with two out 
of phase s-states. In the presence of weak orthorhom- 
bicity, ( and for particular parameter choices) this state 
exhibits the features of a d x 2_ y 2 state, having one gap 
elongated along the a and the other along the b axis of 
the crystal. In this way results compatible with Joseph- 
son experiments on YBaCuO [18] can be obtained in a 
bi-layer s-wave system. Moreover, for this state there will 
be enhancements of the superconducting cross section be- 
low T c along the (tt, 7r) direction, which are rather similar 
to those of the in phase d-wave case discussed above. 

In conclusion, we have applied a theory which has been 
previously extensively compared to normal state mag- 
netic data, in order to ascertain the neutron signatures 
of d-wave pairing. While at low frequencies incommen- 
surate peaks are expected, at higher u> commensurate 
structure will occur associated with a narrow frequency 
peak. The latter feature may have been observed experi- 
mentally by a number of groups. In the present scenario, 
the sharp onset of the frequency peak indicates that there 
is some residual exchange interaction present near opti- 
mal stoichiometry. On the otherhand, the rather broad 
q-space features observed in the normal state for both 
YBa2Cu307 and YBa2Cus06 7 suggests that the mag- 
netism is relatively weak for the highest transition tem- 
perature systems. While this analysis shows consistency 
with d-wave pairing, an alternate orthorhombic bi-layer 
model with two out of phase s- wave bandgaps cannot be 
ruled out. In order to fully verify the origin of the 41 meV 
feature, analogous effects should be searched for in the 
reduced oxygen and LaSrCuO systems. Our calculations 
should provide some guidelines for future experiments on 
these materials. 
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The results of these calculations are in reasonable 
agreement with Ref [5], although the calculated width 
is somewhat wider than observed. The dependence on 
q = (q X) q y ) is similar for both q z = and q z = 7r, but 
the amplitude is only signficiant for latter wave vector. 
This is a consequence of the perpendicular exchange in- 
teraction Jj_ . The mechanism for the formation of the 
strong peak below T c , is as conjectured in Ref. [5], a 
consequence of pair creation with one electron in the + 
lobe and the other electron in the - lobe. However, we 
find that the underlying fermiology of the system also 
plays a role. The separation of the two lobes is near 
(q X) q y ) = (-7T, 7r) . This enhances the weight around 
(q X) q y ) = (-7T, 7r) in the Lindhard function. The further 
inclusion of an antiferromagnetic exchange J(q) creates 
a strong peak at [tt, 7r) . The position of the Fermi 
surface nesting vectors can displace this maximum away 
from (7r, 7r). Thus in LaSrCuO (where there is strong 
nesting close to but not at [tt, 7r) there is a general en- 
hancement of the (q x ,q y ) = (tt, 7r) cross section at twice 
the gap frequency, but the global maxima appear at the 
incommensurate wave-vectors of the normal state neu- 
tron cross section. We estimate that the counterpart of 
these effects should be seen in optimally doped LaSrCuO 
at frequencies of roughly 2A > lAmeV . 

In Figures 2 are plotted the frequency dependences of 
the (7r, 7r) feature (solid line) , with the normal state 
result shown as the dashed line for YBa2Cus07 (a) and 
YBa2Cu306.7 (b). The normal state has a weak max- 
imum in the de- oxygenated case and an even smaller 
feature at optimal stoichiometry, both of which are asso- 
ciated with their respective Van Hove energies. In the ab- 
sence of antiferromagnetic interaction effects, the width 
of the superconducting frequency peak is determined by 
2A at the lower end and the Van Hove energy at the up- 
per. In the presence of J(q), both energies are softened 
by spin fluctuation effects. The gap feature rises more 
sharply the larger is J(q). For the parameters that we 
have chosen from numerous other considerations [6] this 
gives rise to a peak width which is considerably sharper in 
the superconducting state, but which should be close to 
the resolution limited value found experimentally. As the 
oxygen concentration is reduced, these two energy scales 
pull apart and the peak broadens somewhat. Ultimately, 
two distinct frequency features appear, as may have been 
observed by Bourges et al [4] [17]. The presence of two 
bands leads to the two rather weak Van Hove features 
of the figure. It should be stressed that the frequency 
regime over which the normal and superconducting con- 
tributions occur is rather similar. This is a consequence 
of narrow band effects. If the band-widths were taken to 
be broader, the (weak) normal feature will be at higher 
frequencies relative to that of the superconducting state. 
Thus, in this system, also, fermiology plays a role. 
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FIG. 2. Calculated frequency dependence of Imx for (a) 
YBa2Cu3 07 with same parameters as in Fig.l; and (b) pre- 
dictions for YBa2Cu306.7 by assuming A = 15meV. 

Finally in Figure 3, the temperature dependence of 
the peak height at fixed u> = AlmeV is indicated for the 
fully oxygenated sytem. This shows the overall scale of 
the superconducting feature relative to the normal state 
background and its rather dramatic temperature onset. 
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FIG. 3. Calculated temperature dependence of Imx f° r 
YBa2Cu3 07 with same parameters as in Fig.l. 
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, , = X° + Zjxhxli ~ X°i2X?2i)(-Jlq cos g_J - J) 

( J 1 " J(x°n + xh) ~ (J±2iXh + J±i2X°2i) + (x nX°22 - X^iXhKJ 2 ~ J 2 ± ) ' 

Here, Xijihi = 1;2) is the Lindhard function in the layer indices of the bi-layers, and J(q) = — J (cos q x + cosq y ) 
is the antiferromagnetic superexchange within the plane. The antiferromagnetic interlayer coupling Jj_(q) is defined 
as [Jj_i2(q) = —J_i_ e lq * 1 = =^121] which is included in the usual RPA procedure, / is the distance between the closest 
C11O2 layers in YBa2Cu306+ x - The electronic coupling t±_ is embedded in the band structure which enters x° : 

X° = Xn+X? 2 +X2i + X22 (3) 

= 1± f^(xU + x__) + ^7^(x;_ + X_ + ) 

where x° a p{ a jP = +; — ) are the Lindhard functions in the bi-layer bonding and antibonding band indices, and the 
transformation between the layer- and band- representations is: 

Xn=X22= \[X° ++ +X -+X + -+X°- + ] (4) 

x°Y2 = xTi = e ^[x° ++ + X°— ~ X+- ~ X-+] 

Given a satisfactory theory of the normal state data, the implications for a rf-wave superconducting state can be 
readily derived. We previously predicted [6] a low frequency ( < 10 meV) signature associated with a commensurate 
to incommensurate transition below T c . This rf-wave signature has not yet been observed [13,14]. In view of the recent 
experimental focus on higher frequencies, at 41 meV, it is important to extend theoretical calculations to this regime. 

In our numerical analysis we use the same intra-band parameters for YBa2Cus07 and YBa2Cus06 7 as were used 
earlier [6] to address NMR and neutron data. The effects of the superconductivity enter Eq(2) in the usual way [6] 
via the Lindhard function. Figure 1 shows the wave-vector dependence in YBa2Cus07 of the imaginary part of x for 
the d-wave case and in the inset, for the normal metal. In addition we have taken t±_ = 40meV and Jj_ = lOmeV . 
Here the frequency u> = AlmeV is slightly larger than twice the rf-wave gap frequency 2A . Following our solution of 
the bi-layer gap equation, Ref [15], we take the two in- phase, rf-wave gaps in the bi-layer bands to be identical [16]. 




4* 

FIG. 1. Calculated neutron scattering form factor Imx as a function of (q x , q y ) for YBa2Cu3 07. Here it is assumed that the 
superconducting gap A Q = 20meV. 
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We analyze neutron experiments on YBa2Cu3C>6+x at various stoichiometries in 
the superconducting state, within the context of a bi-layer theory which yields good 
agreement with the normal state Cu-NMR and neutron data as a function of to, q 
and P. A d-wave superconducting state exhibits peaks at q = (ir, ir, ir) and sharp 
maxima as a function of w, at twice the gap frequency. Phis behavior may have 
been observed experimentally. Phe counterpart behavior for other choices of order 
parameter symmetry is discussed. 



In view of the growing experimental support for rf-wave 
superconductivity in the YBa2Cu306+ x cuprate family, 
it is becoming increasingly important to characterize the 
magnetic fluctuations in these materials. While a mag- 
netism induced pairing scenario is not uniquely implied 
by this order parameter symmetry, it is entirely consis- 
tent with it. Consequently, neutron experiments with 
their capability of resolving both the q and u> structure 
of the spin fluctuations are of great interest. They pro- 
vide the most direct constraints on the normal state spin 
fluctuations which might drive rf-wave superconductivity 
and in the superconducting state, they contain important 
signatures of the order parameter symmetry. 

In early neutron experiments [1,2], antiferromagnetic 
peaks ( at q = 7T, 7T, 7T ) were observed at reduced oxy- 
gen stoichiometries . As optimal stoichiometry was ap- 
proached, the peak widths were found to become rela- 
tively broad so that even for YBa2Cus06 7 they corre- 
sponded to in plane magnetic correlation lengths which 
were between 1 and 2 lattice spacings. Nevertheless these 
peaks were associated with rather low characteristic en- 
ergy scales ( ~ 20 - 40 meV). Subsequent experiments 
have reported weak [3,1,4] or absent [5] normal state 
magnetic peaks for the fully oxygenated material, and, 
when present, analogously low energy scales. In addition 
a sharp feature at 41 meV has been noted which may [3] 
or may not [5] be associated with weak magnetic scat- 
tering in the normal state. Of particular interest is the 
enhancement of this 41 meV feature below T c . 

Previous work by our group [6,7] has addressed these 
normal state data and explained q, to and temperature T 
dependences as well as trends with variable oxygen stoi- 
chiometry, through a combination of Coulomb renormal- 
ized bandwidth and weak residual exchange interaction 
effects. Similar calculations met with the same degree of 
success for the LaSrCuO family. In contrast to other ap- 
proaches [8,9] it is argued that at optimal stoichiometry, 
these materials are not particularly close to magnetic in- 



stabilities. Incipient localization effects help provide the 
low energy scales. A three band, large U, ( 1/N ) calcu- 
lation yields a dynamical susceptibility 



where, unlike the Hubbard RPA theory, the enhance- 
ment factor contains additional tight binding q-structure, 
J(q) = — J [cosq x a + cosg^a]. The Lindhard function 
X° is given by the usual expression with a three band- 
Coulomb renormalized energy dispersion [6]. As an inde- 
pendent check on the latter, we have confirmed that the 
computed Drude-fitted plasma frequencies [6] are in rea- 
sonable accord with experiment. Our earliest work em- 
phacized [10] the importance of including realistic Fermi 
surface shapes in this analysis, which philosophy has been 
widely adopted by the community. Specifically, for the 
case of YBa2Cu306+ x we find two Fermi surfaces which 
are roughly concentric and which are rotated 45° rela- 
tive to the Fermi surface of the nearest neighbor tight 
binding model. Both of these effects have been observed 
experimentally in photoemission ( ARPES ) experiments 
[11,12]. It should be emphasized that the bandstructure 
does not suggest there are special nesting vectors (e.g. 
(71", 7r)) in YBCO. NMR data, which provide additional 
constraints on models for have also been reasonably 
successfully addressed at the Cu site, although the stan- 
dard hyperfine transfer model for the oxygen relaxation 
appears to be more problematic. 

A bi-layer modulation of the cross section along the 
c-axis is another important feature of the neutron data. 
This was found [7] to be a consequence of antiferromag- 
netic inter- layer correlations Jj_ . While an electronic 
hopping t±_ will produce a c-axis modulation, in order 
for the maximum intensity to occur for all frequencies at 
q z = 7T, this magnetic coupling must be present. To see 
the effects of t±_ and Jj_ we re-write Eq(l) as follows: 
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